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ABSTRACT 

In this paper, we present results from the complete set of cosmic microwave background (CMB) 
radiation temperature anisotropy observations made with the Arcminute Cosmology Bolometer Array 
Receiver (ACBAR) operating at 150 GHz. We include new data from the final 2005 observing season, 
expanding the number of detector-hours by 210% and the sky coverage by 490% over that used for 
the previous ACBAR release. As a result, the band-power uncertainties have been reduced by more 
than a factor of two on angular scales encompassing the third to fifth acoustic peaks as well as the 
damping tail of the CMB power spectrum. The calibration uncertainty has been reduced from 6% 
to 2.1% in temperature through a direct comparison of the CMB anisotropy measured by ACBAR 
with that of the dipole-calibrated WMAP5 experiment. The measured power spectrum is consistent 
with a spatially flat, ACDM cosmological model. We include the effects of weak lensing in the power 
spectrum model computations and find that this significantly improves the fits of the models to the 
combined ACBAR+WMAP5 power spectrum. The preferred strength of the lensing is consistent with 
theoretical expectations. On fine angular scales, there is weak evidence (1-lcr) for excess power above 
the level expected from primary anisotropies. We expect any excess power to be dominated by the 
combination of emission from dusty protogalaxies and the Sunyaev-Zel'dovich effect (SZE). However, 
the excess observed by ACBAR is significantly smaller than the excess power at I > 2000 reported 
by the CBI experiment operating at 30 GHz. Therefore, while it is unlikely that the CBI excess has 
a primordial origin; the combined ACBAR and CBI results are consistent with the source of the CBI 
excess being either the SZE or radio source contamination. 
Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) are among the most powerful and important 
tests of cosmological theory. Measurements of the an- 
gular power spectrum of CMB temperature anisotropies 
on angular sc ales > 10' - corresp onding to multipolcs 
P. < 1000 - (jSoergel et al] 120061) in conjunction with 



other cosmological probes dBurles et al.ll2001b ICole et al.l 
120051: iTegmark et all 120061: iRiess et al.l I2007D have pro- 
duced compelling evidence for the ACDM cosmological 
model. At higher multipolcs, measurements probe the 
Silk damping tail of the power spectrum and provide an 
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independent check of the cosmological model. 

At smaller angular scales, the primary CMB 
anisotropies originating at redshift z = 1100 are exponen- 
tially damped by photon diffusion. This effect, known 
as Silk damping, makes secondary anisotropies - those 
induced along the line of sight at lower redshift - increas- 
ingly important at higher £. At 150 GHz, for example, 
the Sunyaev-Zel'dovich effect (SZE) is expected to be 
brighter than the primary CMB anisotropy at I > 2500. 
The amplitude of the SZE depends sensitively on the am- 
plitude of the matter perturbations, scaling as cr|. Mea- 
surements of the CMB power spectrum with sufficient 
sensitivity on arcminute scales not only extend tests of 
the ACDM model's ability to accurately predict the fea- 
tures in the power spectrum of primary CMB anisotropy, 
but also probe the epoch of cluster formation and provide 
an independent measure of as ■ 

In this paper, we present the complete results of 
observations of CMB temperature anisotropies at 150 
GHz with 5' resolution from the Arcminute Cosmol- 
ogy Bolometer Array Receiver (ACBAR) experiment 
at the South Pole station. Previous measurements 
of the CMB power spe ct rum by ACBAR have been 
presented in iKuo et all ((2004) (hereafter K04) and 
iKuo et all (|2007D (hereafter K07). In addition, the an- 
gular power spectrum on these angular scales has bee n 
meas ured at 30 GHz by C BI dReadhead et alJl2004ari 



VSA (jDickinson et al.ll2004h . and BIMA dDawson et al 



2006), and at 100 and 150 GHz by QuAD (|RAde et al 

Mm- 
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have been consistent with predictions of the primary 
anisotropy based on measurements at larger angular 
scales, with one ex ception. Both CBI ( M ason et al.ll2003t 
iBond et alj I2005D and BIMA (|Dawson et all |2006() ob- 
serve excess power for £ > 2000 at 30 GHz compared 
to the predictions of the ACDM model. This excess can 
be explained by the SZE if erg « 1, but this value is 
in tension with the best-fit WMAP5 value of as ~ 0.8. 
In K07, we found that while the frequency dependence 
of the excess is consistent with the SZE, the ACBAR 
and CBI data could not be used to rule out radio source 
contamination or systematic errors as the source of the 
CBI excess. Careful measurements over a broad range 
of frequencies and angular scales are needed to provide 
a definitive answer. 

Current estimates of the primordial power spectrum 
are consistent with the predictions of slow-roll inflation 
for a nearly scale-invariant spectrum which may also in- 
clude a small running of the spectral index. Sparked 
by the modest evidence for negative running in the 
WMAP first-year data, a number of authors have in- 
vestigated how ex isting data sets limit the allowed infla- 
tionary scenarios (iPeiris et a l. 2003; Muk heriee &: Wangl 
[200l iBridle et al 1 120031: iLeach fc Liddldl2003fh Small 
scale data extend the range over which the primordial 
power spectrum is measured and can potentially yield 
information about the mechanism of inflation. 

This is the third and final ACBAR power spectrum re- 
lease. The first release in K04 analyzed two fields from 
the 2001 and 2002 seasons with a conservative field dif- 
ferencing algorithm. The second ACBAR power spec- 
trum, presented by K07, added two more fields from 
the 2002 season and implemented an improved, undiffcr- 
enced Lead-Main- Trail (no-LMT) analysis of the dataset. 
The results presented here improve on the previous work 
in two ways. First, we include seven additional fields 
observed in the 2005 Austral winter. These fields dou- 
ble the total number of detector hours and substantially 
improve the precision of the band-power estimates. In 
particular, the new fields were selected to dramatically 
expand ACBAR's sky coverage in order to reduce the 
cosmic variance contribution to the uncertainty and to 
improve the multipole resolution on angular scales below 
£ < 1800. This angular range covers the third to fifth 
acoustic peaks, making it especially interesting for con- 
straining cosmological models. Second, we implement 
a new temperature calibration based on a comparison 
of CMB fluctuati ons as measured by ACBAR and the 
WMAP satellite (jHinshaw et all 120081) . This improved 
calibration tightens constraints on cosmological models 
found from the combination of high-^ ACBAR band- 
powers with low-£ results from other experiments. 

This paper is organized as follows. In § [2] we review 
the ACBAR instrument and the CMB observation pro- 
gram. The analysis algorithm is explained in § [3l Section 
§ |4] is an overview of the calibration; the details of cross- 
calibration between WMAP5 and ACBAR are discussed 
in Appendix A. Information on ACBAR's beams can be 
found in § [5] Systematic tests and foreground contami- 
nation are discussed in § [6j We present the band-power 
results in § [3 including a discussion of the scientific inter- 
pretation. The ACBAR band-powers are combined with 
the results of other experiments to place constraints on 
the parameters of cosmological models in § [8] In § [9l we 



summarize the main results of this paper. 

2. THE INSTRUMENT AND OBSERVATIONS 

The ACBAR receiver was designed to take ad- 
vantage of the excellent observing conditions at the 
South Pole to make extremely deep maps of CMB 
anisotropies (|Runvan et al. 1 120031 ). It observes from the 
Viper telescope, a 2.1m off-axis Gregorian with a beam 
size of 5' at 150 GHz. The beams are swept across the sky 
at near-constant elevation by the motion of an actuated 
flat tertiary mirror. 

The receiver contains 16 optically active bolometers 
cooled to 240 mK by a three-stage He 3 -He 3 -He 4 sorption 
refrigerator. The results reported here are derived from 
the 150 GHz detectors: there were 4-150 GHz bolome- 
ters in 2001, 8 in 2002 and 2004, and 16 in 2005. The 
detectors were background limited at 150 GHz with a 
sensitivity of approximately 340 fiKyG. 

In total, ACBAR observed 10 independent CMB fields, 
detailed in Table [1] The power spectrum derived from 
portions of four fields, CMB2/4, CMB5, CMB6, and 
CMB7 was reported in K07. Since then, we have com- 
pleted the analysis of six new fields observed in 2005 as 
well as additional observations of the original four fields. 
Details of the instrument configuration an d performance 
in the 2001 and 2002 seasons are given in lRunvan et al.l 
(2003), while additional details of the CMB observations, 
data reduction procedures, and beam maps can be found 
in K04 and K07. 

3. UN-DIFFERENCED POWER SPECTRUM ANALYSIS 

Following the conventions of the previous data releases, 
the band-powers q are reported in units of [iK 2 , and are 
used to parameterize the power spectrum according to 

£{£+l)C e /2Tr = V e = J2 ( lBXBi , (1) 

B 

where \bi are tophat functions; \bi — 1 for I € B, and 
Xbi = for £ B. The ACBAR observations were car- 
ried out in a lead-main-trail (LMT) pattern. Originally, 
the three fields were differenced according to the formula 
M — [L + T)/2 in order to remove time-dependent chop- 
per synchronous offsets. In K07, this conservative strat- 
egy was shown to be unnecessary and an un-diffcrcnccd 
analysis algorithm was presented. We continued to ob- 
serve in a lead-trail or LMT pattern in 2005 in order to 
produce maps wider than the maximum range (~ 3°) of 
the chopping tertiary mirror. The un-differenced analy- 
sis presented in K07, and used for this paper's analysis, 
is outlined below with any differences in the application 
to the 2005 data set highlighted. 

Let d a be a measurement of the CMB temperature 
at pixel a. The data vector can be represented as the 
sum of the signal, noise and chopper synchronous off- 
sets: d a = s a + n a + o a . For example, although the 
chopping mirror moves the beams at nearly-constant ele- 
vation, the slight residual atmospheric gradient produces 
a chopper synchronous signal o a which is an approxi- 
mately quadratic function of chopper angle. To remove 
these offsets, the data from each chopper sweep are fil- 
tered with the "corrupted mode projection" matrix II to 
produce the cleaned time stream d = lid. 

The II matrix projects out a third to tenth order poly- 
nomial which suppresses large angular scale chopper off- 
sets. The order of the polynomial removed depends on 
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TABLE 1 
CMB Fields 



Field 


RA (dcg) 


dec (deg) 


Area (deg 2 ) 


Year 


# of detectors 


Detector hours 


CMB2(CMB4) 


73.963 


-46.268 


26(17) 


2001(2002) 


4(8) 


2.0k(l.lk) 


CMB5 


43.371 


-54.836 


28 


2002(2005) 


8(16) 


13.2k(10.6k) 


CMB6 


32.693 


-50.983 


23 


2002 


8 


2.8k 


CMB7(ext*) 


338.805 


-48.600 


28(107) 


2002(2005) 


8(16) 


3.4k(12.4k) 


CMB8 


82.297 


-46.598 


61 


2005 


16 


17.1k 


CMB9* 


359.818 


-53.135 


93 


2005 


16 


4.0k 


CMB10* 


19.544 


-53.171 


93 


2005 


16 


3.6k 


CMB11* 


339.910 


-64.178 


91 


2005 


16 


4.9k 


CMB12* 


21.849 


-64.197 


92 


2005 


16 


2.7k 


CMB13* 


43.732 


-59.871 


78 


2005 


16 


7.6k 



Note. — The central coordinates and size of each CMB field observed by ACBAR. The sixth column 
gives the number of 150 GHz detectors. The last column gives the detector integration time for each 
field after cuts. The detector sensitivity was comparable (within ~10%) between 2002 and 2005. The six 
largest fields (marked with a *) are used in the calibration to WMAP. Note that the 2005 observations 
extended the declination range of the CMB7 field, leading to the combined field CMB7ext. CMB2(CMB4) 
and CMB8 also partially overlap, but are analyzed separately for computational reasons. Approximately 
1/4 of the CMB2(CMB4) scans have been discarded to eliminate the overlapping coverage. The listed 
numbers reflect this loss. 
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Fig. 1. — The ACBAR fields overlaid on the IRAS dust map. The position of each field is plotted and labeled with the field name. The 
color coding indicates the year in which the observations occurred: red = 2001, orange = 2002, and yellow = 2005. The bulk of the 2005 
season was targeted at large, co mparatively shallow fields, increasing the total sky coverage by a factor of six. The fields are plotted on top 
of the 100 fim IRAS dust map l|Schlegel et aljri998ft . Each field lies within the "Southern hole", a region of low dust emission visible from 
the South Pole. The CMB8 field (lower right corner) was targeted at the deep region of the B03 experiment as an alternative calibration 
path to the WMAP cross-calibration used for the results presented here. 

the amplitude of atmosphere-induced cross-channel cor- 
relations. As described in K07, small angular scale offsets 
can be be removed by subtracting an "average" chopper 
function. For 2002 data, we removed a chopper syn- 
chronous offset from each data strip where the ampli- 
tude of the offset at each sample in the strip is free to 
vary quadratically with elevation in the map. The large 
fields observed in 2005 have up to four times the dec 
range of the fields observed in 200f and 2002 (~ 10° vs. 
~ 2.5°). For 2005 data, we allow the offset to vary from 
a third to fifth order polynomial depending on the extent 



of the map in declination. A zeroth order polynomial in 
elevation removes the average chopper function and the 
higher order terms effectively act as a high-pass filter 
on changes in the offset as a function of time or eleva- 
tion. This anisotropic filtering removes offset-corrupted 
modes while preserving most of the uncorrupted modes 
for the power spectrum analysis. The loss of information 
at high-£ is small; the removed modes account for only a 
few percent of the total degrees of freedom of the data. 

The corrupted mode projection matrix LT can be rep- 
resented as the product of two matrices, II = n 2 ITi. 
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The operator 111 is the original II matrix referenced in 
K04 which adaptively removes polynomial modes in RA. 
The additional operator II2 removes modes in dec in- 
dependently for each of the lead, main, and trail fields 
and can be further decomposed into the product II2 = 
j\P ol VY\LPF rp ne p era ^ or n^ ol y performs the aforemen- 
tioned polynomial projection in dec to remove small- 
scale chopper offsets. The second operator H2 PF im- 
poses a low-pass filter (LPF) I < 3200 on each dec strip. 
The dec strips are perpendicular to the scan direction; 
the timestreams have always had a LPF applied in the 
scan direction. The pixelation used when estimating the 
power spectrum is too large to resolve all of the noise 
power (at I up to 10,800), causing out-of-band noise to 
be aliased into the signal band (I < 3000) if a LPF is not 
applied. Eliminating this high-frequency noise reduces 
the contribution of instrumental noise to the reported 
band-powers. 

Using the pointing model, the cleaned timestreams are 
coadded to create a map: 

A = Ld. 

The noise covariance matrix of the map can be repre- 
sented as 

C N = L(nn')L*. 

where n is the timestream noise. The noise matrix is 
diagonalized as part of applying a high signal-to-noise 
transformation to the data. Eliminating modes with 
insignificant information content reduces the computa- 
tional requirements of later steps in the analysis. 

In order to apply the iterative quadratic band-power 
estimator, we need to know the partial derivative 
of the theory covariance matrix with respect to each 
band-power qs- The theory matrix can be calculated 
by considering the effects of the filtering on the raw sky 
signal. The signal timestream s a is the convolution of 
the true temperature map T(r) with the instrumental 
beam function B a (r) 

s a = J d 2 r1(r)B a (r). 

The signal component of the coadded map will be 
A si s = Ls or 

AT 9 = J d 2 rF l (r)T(r), 

where we have defined the pixel-beam filter function Fi 
fi(r) = 5^L fa B a (r). 

a 

The theory covariance matrix can be calculated in the 
flat sky case to be 

C T{ij} ee (A t A 3 r 9 =// d 2 rd 2 r'F 4 (r)^(r')(2:(r)T(r')) 




.il-(r-r') 



= /(0 Q -^ (1) ^ (1) ' ^ 



where is the Fourier transform of the pixel-beam 

filter function. The partial derivative of the theory ma- 
trix can be calculated in a straightforward manner from 
equations [1] and [2j 

This algorithm does not require the instrument beams 
to remain constant. The actu al ACBAR beam siz es vary 
slightly with chopper angle (|Runvan et al.l [20031. The 
measured beam variations can be fit to a semi-analytic 
function as described in K04 to create a more accurate 
representation of the true beam shape across the map. 
We use the corrected beam sizes when removing point 
sources. In K04 and K07, we found that the differences 
in the power spectra from using the map-averaged beam 
or exact beam for each pixel were negligible. For the 
band-powers reported in Table [3l an averaged beam is 
used for the entire map. 

As in K07, we calculate the full two dimensional noise 
correlation matrix directly from the time stream data 
without using Fourier transforms. All the numerical 
calculations are performed on the National Energy Re- 
search Scientific Computing Center (NERSC) IBM SP 
RS/6000. The evaluation of the Fourier transform of 
Fi (r) is the most computationally expensive step of this 
analysis. We use an iterative quadrati c estimator to fin d 
the maximum likelihood band-powers (|Bond et al.|[T9 98V 
The resulting band-powers are presented in Table [3] and 
Figure 

4. CALIBRATION 

We derive the absolute calibration of ACBAR by di- 
rectly comparing the 2005 ACBAR maps to the WMAP 5 
V and W-band temperature maps (jHinshaw et al.| [2008). 
We pass the WMAP5 maps through a simulated ver- 
sion of the ACBAR pipeline to ensure equivalent filter- 
ing and cross-spectra are calculated for each field. The 
ratios of the cross-spectra are used to measure the rel- 
ative calibration after being corrected for the respective 
instrumental beam functions. We had initially applied 
this calibration scheme to the WMAP3 maps. Transi- 
tioning to the WMAP5 dataset lowered the calibration 
by 1.4% in CMB temperature units and slightly reduced 
the overall uncertainty. The ACBAR band-powers are 
unchanged except for this calibration factor. Results for 
ACBAR's six largest fields (approximately 600 deg 2 in 
area) are combined to achieve a calibration accuracy of 
1.97% for the 2005 data. 

The 2005 calibration is transfered to 2001 and 2002 
through a comparison of power spectra for overlapping 
regions observed by ACBAR in each year. The CMB5 
field is used to extend the calibration of the 2005 sea- 
son to the 2002 data. The CMB5 calibration is carried 
to other fields observed in 2002 by daily observations of 
the flux of RCW38. The calibration of the CMB4 field 
(observed in 2002) then is transfered to the 70% overlap- 
ping CMB2 field (observed in 2001). According to the 
new calibration, results from the RCW38-based calibra- 
tion used for the 2002 data in K07 need to be multipled 
by 0.959 ± 0.032. Including the year-to-year calibration 
uncertainty, the final calibration has an uncertainty of 
2.05% in CMB temperature units (4.1% in power). Addi- 
tional details of this procedure are discussed in Appendix 

El 

5. BEAM DETERMINATION 
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The beams are well-described by a symmetric Gaus- 
sian, with their main-lobe FWHM determined to 2.6% by 
continuous measurements of the images of bright quasars 
located in the CMB fields. The beam sidelobes were mea- 
sured to the level of 30 dB with observations of Venus 
made in 2002. Venus is extremely bright at millimeter 
wavelengths and with a diameter of < 1', is much smaller 
than ACBAR's beam size. However, there are extended 
periods during which ACBAR was unable to observe 
Venus. ACBAR observed RCW38, a bright HII region in 
the galactic plane, every day. We compare deep, coadded 
observations of RCW38 to constrain the temporal vari- 
ability of the beam sidelobes when Venus was unavail- 
able. The complex structure surrounding RCW38 makes 
it difficult to directly recover the beam shape B(r). In- 
stead, we monitor ratios of the beam-smoothed RCW38 
maps J d 2 rS RCW38 (r)B(r). Any observed differences in 
the maps would indicate a change to the instrumental 
beam function as the morphology of RCW38's emission 
S RCW3S is expected to be constant. We set an upper 
limit on the possible temporal variations in the map and 
use this to constrain temporal variations in the beam 
function. The estimated band-power uncertainty from 
the beam function is comparable to the overall calibra- 
tion uncertainty and is plotted in Figure [2] 

6. SYSTEMATIC UNCERTAINTIES AND FOREGROUNDS 

6.1. Jackknife Tests 

We performed a scries of tests to search for and con- 
strain potential systematic errors in the power spectrum 
results. As described in K04, the data can be divided 
into two halves based on whether the chopping mirror is 
moving to the left or right. The "left minus right" jack- 
knife is a sensitive test for errors in the transfer function 
correction, microphonic vibrations excited by the chop- 
per motion, or the effects of wind direction. Maps with 
bright sources such as RCW38 can provi de particularly 
sensit ive tests of the transfer function (see lRunyan et all 
(2003) for a description of ACBAR's transfer functions). 
Similarly, the data can be split based on the time that 
the observation occurred. A non-zero signal could be pro- 
duced in the "first half minus second half" jackknife by 
variation in the calibration, pointing, beam and sidelobe, 
or any other time dependent variations in the instrument. 
In addition, the band-powers of each jackknife constrain 
the mis-estimation of noise during that period. 

We applied the left-right jackknife to the 2005 CMB 
data and found the band-powers were inconsistent with 
zero at 2.5c at high-^ (£ > 2100). We reran a set of left- 
right jacknives dropping individual channels, and found 
that two channels stood out. With both channels ex- 
cluded, the discrepancy in the left-right jackknife band- 
powers disappeared. We were unable to find evidence 
for unusual microphonic lines or transfer functions in the 
two problematic channels, but hypothesize that these two 
channels have subtle microphonic response in the signal 
bandwidth that are detectable only in a deep integra- 
tion. Both channels are excluded from the 2005 data for 
all results reported in this paper. 

We apply the first-second half jackknife test to the 
joint CMB power spectrum with the exclusion of the 
bad channels from the 2005 data. The power spectrum 
of the chronologically differenced maps is compared to 
the band-powers of a set of Monte Carlo realizations of 



simulated difference maps in order to account for a num- 
ber of effects such as the small filtering differences due 
to different scan patterns and the temporal uncertainty 
in the beam sidelobes (see We find that the jack- 
knife band-powers are consistent with the predictions of 
the Monte Carlo above £ = 400. There is a 4cr resid- 
ual of ~15 /iK 2 in the first bin. Because the combined 
statistical and cosmic variance uncertainty in this bin is 
a factor of six larger, we assume that the band-power 
estimate will not be significantly biased. 

We also perform the left-right jackknife on the joint 
CMB power spectrum found from the complete data set. 
The results are consistent with zero for £ > 900. Sta- 
tistically, the probability to exceed the measured x 2 for 
£ > 900 is 15%. The results are inconsistent with zero 
at a very low (~4 /J.K 2 ) level (Fig. [3]) on larger angu- 
lar scales. This residual could be due to a small noise 
mis-estimate at low-^, possibly caused by neglected at- 
mospheric correlations. The jackknife failure of ~4 /iK 2 
is much smaller than the band-power uncertainties (90 - 
300 ^K 2 ) in these £-h\as which are dominated by cosmic 
variance. The first-second half jackknife is insensitive to 
discrepancies of this magnitude due to the greater uncer- 
tainties introduced by small pointing and filtering differ- 
ences. We conclude that the complete ACBAR data set 
shows no significant residuals in the jackknife tests and 
we expect no significant systematic contamination of the 
resulting power spectrum. 



6.2. Foregrounds 

The potential contribution of foreground emission 
must be considered in the interpretation of CMB tem- 
perature anisotropics. There are three foregrounds with 
significant emission at 150 GHz on the relevant angu- 
lar scales: galactic dust, extragalactic radio sources, and 
dusty proto-galaxies. As an effectively single- frequency 
instrument, ACBAR depends on data from other ex- 
periments to construct and constrain foreground mod- 
els. We use the methodology described in K04 to re- 
move templates for radio sources and dust emission from 
the CMB maps without making assumptions about their 
flux. The contribution from dusty protogalaxies is less 
certain; however, we do not expect the combined resid- 
ual foreground emission to significantly impact the power 
spectrum for £ < 2400. 

We remove modes from the CMB maps correspond- 
ing to radio sou rces in the 4.85 GH z Parkes-MIT-NRAO 
(PMN) survey (jWright et al.l Il99l . Of the 1601 PMN 
sources with a flux greater than 40 mJy that lie in the 
ACBAR fields, we detect 37 sources including the guid- 
ing quasars at greater than 3er with the application of 
an optimal matched filter. There are less than 2.2 false 
detections expected with this detection threshold. The 
measurement errors are estimated through sampling the 
distribution of pixels in a set of 100 Monte Carlo real- 
izations of the CMB+noise for each field. Table [2] lists 
the parameters of the detected PMN sources. Except 
for the few bright sources detected at 150 GHz, remov- 
ing the PMN point sources does not significantly affect 
the band-powers. 

It is possible that faint radio sources, undetected at 
150 GHz, could contribute to the observed band-powers. 
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Fig. 2. — ACBAR Beam Uncertainty and Beam Function. Solid line & left axis: The la envelope for uncertainty in the ACBAR beam 
function Bg . The increasing uncertainty above I = 1000 reflects the 2.6% uncertainty in the fitted Gaussian FWHMs. The behavior at 
t < 1000 is a combination of the uncertainty in the measured sidelobes and the calibration method 'pinning' the transfer function for 
t £ [256, 512], Dashed line & right axis: The measured ACBAR beam function. 
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Fig. 3. — Systematic tests performed on the ACBAR data. Top: Power spectrum (red triangle) for differenced maps from the first 
half of the season and second half of the season for each field, compared to the results of Monte Carlo simulations (error bars). Middle: 
Power spectrum (blue star) derived from difference maps of the left- and right-going chopper sweeps for all ten fields. Bottom: The 
undifferenced band-powers from Table \3\ (black diamonds) compared to both jackknife power spectra: the left-right jackknife (blue star) 
and first half-second half jackknife (red triangle). 
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We parameterize the contribution as 

2 



T> SIC — a 



2600 



(3) 



which is appropriate for unclustered point sources. We 
compare the 150 GHz ACBAR point source n umber 
counts to the model in lWhite &: Maiumdarl (|2004f ) based 
on WMAP Q-band data, 

r/.V _ 80 deg- 2 ' '"' x " J ; 
dS v 



1 mJy \ 1 mjy / 

to constrain the residual power contribution at 150 GHz. 
We can use this model to estimate the residual band- 
power contribution due to sources too faint to be included 
in the PMN catalog. 

Following the convention in that work, the spectral de- 
pendence of the fluxes is parameterized as S v cx v^ 3 . The 
number of PMN sources detected at 150 GHz in a log- 
arithmic flux bin, n# , are compared to the predicted 
number counts from the same population of sources with 

a given /3, + rig olse . Here, rijP is the modeled num- 
ber counts and ng mse is the expected number of false 
detections due to ACBAR's measurement uncertainty. 
The number counts are assumed to follow a Poisson dis- 
tribution. Sources with estimated measurement errors 
greater than 140 mJy in the ACBAR maps are cut to 
reduce the n 7 jf lse term. The modeled number counts 

are scaled by (N to t ~ N cut )/N t ot to compensate. All 
other sources with measured amplitudes greater than 350 
mJy at 150 GHz are included in the calculation with- 
out consideration of the signal-to-noise. If the sources 
found by ACBAR at 150 GHz are the same population 
found by WMAP, this implies a uniform spectral index 
of/3 = 0.14±0.15. However, this small sample selected 
for high flux at 150 GHz is heavily biased toward sources 
with flat or rising spectra. We increase ACBAR's sensi- 
tivity to dimmer sources by binning all sources within a 
given PMN flux range, and look at the ratio of the av- 
erage flux at 150 GHz to the average flux at 4.85 GHz 
within each bin. We find the ratio (Siso/Skss) increases 
with PMN flux from 0.07 for sources below 400 mJy to 
0.41 for sources above 1600 mJy at 4.85 GHz. This im- 
plies that the sources in the PMN catalog have a flux- 
dependent spectral index where dimmer objects typically 
have a more steeply falling spectrum. The band-power 
contribution of the low-flux sources depends sensitively 
on the extrapolation of the 40 mJy flux cutoff in the 
4.85 GHz PMN catalog to 150 GHz. Based on the ob- 
served flux ratios for PMN sources with S4.85 < 400 mJy, 
we conservatively assume S150/S4.&5 —0.1 for a flux cut- 
off at 4mJy at 150 GHz. This flux ratio corresponds to a 
spectral index of (3 = —0.67, well below f3 = 0.14 ± 0.15 
estimated from the ACBAR detected sources and the 
WMAP Q-band source model. Estimating the residual 
radio source band-power contribution at 150 GHz with a 
flux cutoff of 4mJy gives g™ dio = 2.2. At this level, the 
residual contribution from radio sources will be negligible 
in the ACBAR data. 

The ACBAR fields are positioned in the "South- 
ern Hole," a region o f exceptionally low Gal actic dust 
emission (Figure [J). iFinkbeiner etafl l(l999) (FDS99) 
constructed a multi-component dust model that pre- 
dicts thermal emission at CMB frequencies from the 



combined observations of IRAS, COBE/DIRBE, and 
COBE/FIRAS. Taking into account the ACBAR filter- 
ing, the FDS99 model 13 predicts a RMS dust signal at 
the few \xK level primarily on large angular scales. The 
ACBAR maps can be decomposed as the sum of the 
CMB and dust signals Tqmb + £,Tp£>s- The dust am- 
plitude parameter £ is predicted to equal unity by the 
FSD99 model. The ACBAR maps are cross-correlated 
with the dust templates Tp^s to calculate the ampli- 
tude in each field. The errors are estimated by applying 
the same procedure to 100 CMB+noise map realizations 
for each field. The uncertainty in £ is dominated by CMB 
fluctuations. The best-fit amplitude from combining all 
the fields is£ = 0.1±0.5. The estimated amplitudes of 
the individual fields are shown in Figure [4] The reduced 
X 2 of the measured amplitudes £s of the eight fields an- 
alyzed is 0.75 for the no-dust assumption of (£) = and 
increases to x 2 — 1-12 for the FDS99 model amplitude 
of (0 = 1- Therefore, the ACBAR data slightly favor 
a lower amplitude than predicted by the FDS99 model. 
The dust signal is not detectable in any of the ACBAR 
fields, and removing the dust template has a negligible 
impact on the measured power spectrum. 

Dusty IR galaxies are the third and least con- 
strained foreground in the ACBAR fields. This pop- 
ulation of high-redshift, star- forming galaxies has been 
studied by se veral experiments at higher frequencies 
(ICoppin et"aT]|2006t lLaurent et all 120051 : iMalonev etHI 
l2005HGreve et al.ll2004( l. However, as discussed in K07, 
extrapolating the expected signal to 150 GHz remains 
highly uncertain, and there remain significant uncertain- 
ties in the number counts ^ and spatial clustering of 
the sources. The frequency dependence can be empir- 
ically determined by comparing the measured number 
counts in overlapping fields observed at different frequen- 
cies. This comparison has been done with MAMBO (1.2 
mm) and SCUBA (850 ^m), leading, to a spectral de- 
pendence of S v oc v 65 (|Greve et al.ll2004h . A second 
method of estimating the index used nearb y galaxy data 
to obtain S u ~ ^2.6fo.3 (|Knox et al.ll2004fl . The uncer- 
tainty in the spectral dependence significantly affects the 
extrapolation of the flux of dusty galaxies to 150 GHz. 
We use estimates of the source numbe r counts from 
the SHADES survey (ICoppin et al.ll2006D an d Bolocam 
Lockman Hole S urvery (IMalonev et al.ll2005ft . We apply 
the formulas in IScott fc White! ( 19990 to estimate the 
expected power spectrum for the source number counts, 
ignoring the clustering terms. In this limit, T>g will have 
the form in eq. [3l Scaling the results to 150 GHz with 
the MAMBO/SCUBA prescription of S v cx u 2S5 leads 
to an estimated contribution of qf^ c sty — 17 — 29. This 
range reflects the differences between the measured num- 
ber counts, but does not include the uncertainty in the 
spectral dependence of the fluxes. Combining the me- 
dian index with earlier SCUBA d ata fit by two power 
laws in S from (|Borvs et al.ll2003f) . we find an excess of 
22 [1K 2 at I = 2600, within that range. This level is only 
a factor of two smaller than the instrumental noise of 
ACBAR and might influence the interpretation of high-^ 
band-powers. We tentatively assume that contamination 
from dusty proto-galaxies does not significantly effect the 
resulting power spectrum. The implications of relaxing 

13 We use the default model 8 of FDS99. 
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this assumption for cosmological parameter estimation 
are explored in § 18.31 

7. RESULTS AND DISCUSSIONS 
7.1. Power Spectrum 

The power spectrum presented in Figure [5] is produced 
by the application of the analysis algorithm outlined in 
§[3] to ACBAR data from the 2001, 2002 and 2005 aus- 
tral winters. The resulting power spectrum is compared 
to the WMAP5 and B03 spectra in Figure El The zero- 
curvature, ACDM "ACBAR+WMAP5" best fit model 
is shown in each figure for reference. The decorrelated 
band-powers are tabulated in Table [3] Our choic e of th e 
decorrelation transformations follows iTegmarkl (|1997f ). 
The band-powers can be com pared to a t heoretical model 
using the window functions (Knoxl [l~999l ). As in K04, we 
sample the likelihood function £(A) = -J= e -( A<c ~ lA )/ 2 

near the maximu m and fit the res ults with offset log- 
normal functions (jBond et al.ll2000T ). The fit parameters 
cr, x are listed in Table[3]as well. The band-powers, likeli- 
hood fit parameters, and window functions are available 
for download from the ACBAR website 14 . 

The ACBAR data extend the measurement of the tem- 
perature anisotropics well into the damping tail with 
S/N > 5 for I < 2300. The fourth and fifth acous- 
tic peaks are detected for the first time in the ACBAR 
band-powers, providing a dditional support for the coher- 
ent origin of anisotropy (|Albrecht et al.| [l996). The po- 
sition of the third acoustic peak is consi stent with previ- 
ous de tections of the feature by CBI (|Rcadhca d et al.l 
l2004ah. B03 Hones et al.1 12006). ACBAR (K07), and 
QUaD (|P. Ade et al.ll2007f) . The ACBAR band-powers 
are in excellent agreement with the cosmological mod- 
els constrained by observations on larger angular scales. 
The probability to exceed the reduced x 2 between the 
ACBAR band-powers and the WMAP3 only best-fit 
ACDM model is 17%. This probability increases to 62% 
with the WMAP5 only best-fit model. This serves as 
both a powerful confirmation of our basic cosmological 
model and an indication of the quality of the ACBAR 
data set. 

7.2. Anisotropics at £ > 2000 

Several theoretical ca lculations (jCoorav eTaTl [20001 : 
iKomatsu &: Seliakl 12002) and hydrodyna mical simula- 
tions (|Bond et all l2005t 1 White et aLll2002ft suggest that 
the thermal Sunyaev-Zel'dovich effect power spectrum 
will exceed that of the primary CMB temperature 
anisotropies for I > 2500 at 150 GHz. The amplitude 
of the SZE power spectrum is closely related to the am- 
plitude of matter perturbations which is commonly pa- 
rameterized as < 7s; the SZE power spectrum is expected 
to scale as erj ([Zhang et al.ll2002t ). To a lesser extent, 
the level of the SZE will also depend on details of clus- 
ter gas physics and thermal history. The non-relativistic 
thermal SZE (AT$z) has a unique frequency dependence 



AT. 



sz 



where x = 



Tcmb 

hv 



y 



(4) 



z^/56.8GHz. The variable y is the 



Compton parameter and is proportional to the integrated 

14 http:/ /cosmology.berkeley.edu/group/swlh/acbar/index.html 



electron pressure along the line of sight. The CBI ex- 
tended mosaic observations (Rcadhc ad et al.ll2004al ) de- 
tected more power above £ — 2000 than is expected from 
primary CMB anisotropies. This excess could be the 
first d etection of the SZE power spectrum (Maso rTet al.l 
120031 iReadhead et~aT1 l2004al : iBond et all l2005h How- 
ever, there are alternative explanations for the ob- 
served power ranging from an unresolved population 
of low-f lux radio sources to non-standard inflationary 
models (iCoorav fc MelclfiOTrl [2u0l i Griffiths et al.l l2Q03l: 
ISubramanian et al.ll2003h that produce higher than ex- 
pected CMB anisotropy power at small angular scales. 
The frequency dependence of the excess power can be ex- 
ploited to help discriminate between the SZE and other 
potential explanations. 

The ACBAR band-powers reported in this paper are 
slightly larger at £ > 2000 than expected for the 
"ACBAR+WMAP5" best fit model. We subtract the 
predicted band-powers at £ > 1950 from the measured 
band-powers in Tableland find an excess of 22 ±20 /j,K 2 
in a flat band-power from 1950 < £ < 3000. This es- 
timate ignores the band-power contribution from dusty 
proto-galaxies which is expected to be comparable (see 
gOJ). The ACBAR band-powers at 150 GHz can be 
used to place constraints on frequency spectrum of the 
larger CBI excess measured at 30 GHz. We parameter- 
ize the excess power for £ > 1950 at the two frequencies 
as P30 = aPi5o and sample the likelihood surface for 
a e [0,10] and P 15Q S [0,300] [iK 2 . The beam uncer- 
tainty and the calibration error for both experiments is 
taken into account by Monte Carlo techniques. The like- 
lihood function is averaged over 1000 realizations under 
the assumption that each of the errors has a normal dis- 
tribution. The resulting likelihood function for a (after 
P150 is marginalized) is shown in Figure [7] From the 
ACBAR and CBI frequency bands, we expect a = 4.3 
for power originating from the SZE. If the excess is due 
to primary CMB anisotropies, we expect a = 1. We con- 
clude that it is more than 5 times as likely that the excess 
seen by CBI and ACBAR is caused by the thermal SZE 
than a primordial source. We expect the contribution of 
radio sources to CMB power to be at least a factor of 
ten higher at 30 GHz than at 150 GHz (a > 10). Be- 
cause of the relatively weak detection of excess power 
by ACBAR, flat spectrum radio sources are determined 
to be ~10% more likely than the SZE to be the source 
of the excess. The lower level of excess power seen by 
ACBAR argues against the the CBI excess having a pri- 
mordial origin, but is consistent with either the SZE or 
radio source foregrounds. When we include the expected 
contribution of dusty protogalaxies to the ACBAR ex- 
cess band-power, the likelihood of the CBI excess being 
due to radio sources increases with respect to thermal 
SZE. 

8. COSMOLOGICAL PARAMETERS 

8.1. Cosmological Parameters and their "Prior" 
Measures 

In this section, we estimate cosmological parameters 
for a minimal inflation-based, spatially-flat, tilted, gravi- 
tationally lensed, ACDM model characterized by six pa- 
rameters, and then investigate models with additional 
parameters to test extensions of the theory. For our base 
model, the six parameters are: the physical density of 
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TABLE 2 
Millimeter Bright PMN Sources 



Source Name/Position 


Field 


Sd «r frnJv) 


5"i fmJv) 


OL 1 Kd 1 A QK 
lOU/ 4.00 


PMN J0455-4616* 


CMB2 


1653 


2898 ± 60 


0.16 


PMN J0439-4522 


CMB2 


634 


383 ± 73 


-0.15 


PMN J0451-4653 


CMB2 


541 


360 ± 58 


-0.12 


PMN J0253-5441* 


CMB5 


1193 


1260 ± 62 


0.02 


PMN J0223-5347 


CMB5 


397 


176 ± 27 


-0.24 


PMN J0229-5403 


CMB5 


242 


147 ± 17 


-0.14 


PMN J0210-5101* 


CMB6 


3198 


1268 ± 86 


-0.27 


PMN J2207-5346* 


CMB7ext 


1410 


381 ± 68 


-0.38 


PMN J2235-4835* 


CMB7ext 


1104 


1509 ± 75 


0.09 


PMN J2239-5701* 


CMB7ext 


1063 


501 ± 68 


-0.22 


PMN J2246-5607 


CMB7ext 


618 


386 ± 51 


-0.14 


PMN J2309-5703 


CMB7ext 


56 


257 ± 76 


0.44 


PMN J0519-4546* 


CMB8 


15827 


1375 ± 101 


-0.71 


PMN J0519-4546* 


CMB8 


14551 


1148 ± 86 


-0.74 


PMN J0538-4405* 


CMB8 


4805 


7114 ±87 


0.11 


PMN J0515-4556* 


CMB8 


990 


671 ± 96 


-0.11 


PMN J0526-4830 


CMB8 


425 


82 ±25 


-0.48 


PMN J0525-4318 


CMB8 


217 


99 ±25 


-0.23 


PMN J0531-4827 


CMB8 


142 


96 ±25 


-0.11 


PMN J2357-5311* 


CMB9 


1782 


347 ± 49 


-0.48 


PMN J2336-5236 


CMB9 


1588 


233 ± 57 


-0.56 


PMN J2334-5251 


CMB9 


557 


432 ± 57 


-0.07 


PMN J0018-4929 


CMB9 


142 


178 ± 57 


0.07 


PMN J0026-5244 


CMB9 


40 


192 ± 62 


0.46 


PMN J0050-5738* 


CMB10 


1338 


773 ± 108 


-0.16 


PMN J0058-5659* 


CMB10 


739 


514 ±61 


-0.11 


PMN J0133-5159* 


CMB10 


672 


248 ± 72 


-0.29 


PMN J0124-5113* 


CMB10 


308 


335 ± 48 


0.02 


PMN J2208-6404 


CMB11 


53 


136 ± 44 


0.27 


PMN J0103-6438 


CMB12 


395 


268 ± 65 


-0.11 


PMN J0144-6421 


CMB12 


152 


184 ± 59 


0.06 


PMN J0303-6211* 


CMB13 


1862 


423 ± 63 


-0.43 


PMN J0309-6058* 


CMB13 


1103 


596 ± 80 


-0.18 


PMN J0251-6000 


CMB13 


433 


189 ± 34 


-0.24 


PMN J0236-6136 


CMB13 


406 


365 ± 33 


-0.03 


PMN J0257-6112 


CMB13 


178 


104 ± 33 


-0.16 


PMN J0231-6036 


CMB13 


174 


105 ± 33 


-0.15 



Note. — These sources from the PMN 4.85 GHz catalog are detected at > 3.0<r sig- 
nificance with ACBA R, corresponding to a false detection rate of < 2.2. The fluxes at 
4.85 GHz (S 4 .85, from I Wright et al.1 lfl99l )) and 150 GHz (ftso, measured by ACBAR) are 
given. For ACBAR, the flux conversion factor is 1 \iKcmb = 0.9 mJy. The spectral index 
a is defined as S v oc v a . The flux of some of these sources varied by up to 50% between 
years; this variability is not reflected in the estimated errors. The central guiding quasars 
(one in each of the 5 deeper fields) are marked with asterisks (*). These sources, as well 
as all other PMN sources > 40mJy, are projected out from the data using the meth- 
ods described by K04 and do not contribute to the power spectrum results. The brightest 
sources are marked with circles (°) and are removed from the maps in a beam- independent 
method. Note that PMN J0519-4546a/b are within one beam width of each other and are 
not separately resolved by ACBAR. As a result, the listed a for PMN J0519-4546a/b is 
estimated from the sum of the fluxes at 4.85 GHz and the mean of the fluxes at 150 GHz. 



baryonic and dark matter, fl^h 2 and £l c h 2 ; a uniform 
spectral index n s and amplitude In A s of the primordial 
power spectrum, the optical depth to last scattering, t; 
and 0, the ratio of the sound horizon at last scattering 
to the angular diameter distance. The primordial co- 
moving scalar curvature power spectrum is expressed as 
V B (k) = A s {k/k i ,) <Jls ^ 1 \ where the normalization (pivot- 
point) wavenumber is chosen to be k* = 0.05 Mpc -1 . 
The parameter maps angles observed at our location 
to comoving spatial scales at recombination; changing 
shifts the entire acoustic peak/ valley and damping pat- 
tern of the CMB power spectra. Additional parameters 
are derived from this basic set. These include: the energy 
density of a cosmological constant in units of the critical 
density, 17a; the age of the universe; the energy density 
of non-relativistic matter, f2 m ; the rms (linear) matter 



fluctuation level in 8/1 Mpc spheres, erg; the redshift to 
reionization, z re ; and the value of the present day Hubble 
constant, Ho, in units of km s -1 Mpc -1 . 

Single-field models of inflation predict the existence 
of a gravitational wave background characterized by a 
primordial power law Vt ~ k nt . We characterize the 
strength by the tensor-to-scalar ratio r = Vt/V s evalu- 
ated at a pivot point 0.002 Mpc -1 . We relate the tilt 
to r using the approximate consistency relation nt ~ 
— r/8/(l — r/16). (We find little difference in the pa- 
rameters if we just fix n t to be zero, as has often been 
assumed when r is included, but using this relation is 
superior since it is motivated by inflation physics.) 

A small running of the spectral index is also expected 
in slow-roll inflation and we test for this by extending 
the basic ACDM power law model to include a scale de- 
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Fig. 4. — Dust emission is not detected in the ACBAR fields. Parameterizing the dust signal as Tcmb + £,TfbS> a suite of Monte Carlo 
realizations of maps of CMB and noise is used to estimate the uncertainty in §. We find the upper limits in each field to be consistent with 
the FDS99 model (£ = 1), but the data somewhat favor a lower dust amplitude. The reduced x 2 of the measured amplitudes £s is 0.75 
under the assumption that (£) = (the dashed line). The reduced \ 2 f° r the FDS99 model with (£) = 1 is 1.12 (the dotted line). 
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Fig. 5. — The decorrelated ACBAR band-powers for the full data set. The lcr error bars are derived from the offsct-lognormal fits to 
the likelihood function. The band-powers are in excellent agreement with a ACDM model. The damping of the anisotropies is clearly seen 
with a S/N > 4 out to £ = 2500. The third acoustic peak (at £ ~ 800), fourth acoustic peak (at i ~ 1100), and fifth acoustic peak (at 
I ~ 1400) are visible. The plotted lines are the best fits to the ACBAR and WMAP5 band-powers for a spatially flat, ACDM universe 
with no SZE contribution. A lensed (red) and unlensed (blue) model spectrum is shown for a fixed parameter set; the lensed spectrum is a 
significantly better fit to the ACBAR data. 



pendence of the scalar spectral tilt, dn s /d\ii{k). 



We have also added non-zero curvature flk to our ba- 
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TABLE 3 
Joint Likelihood Band-powers 



i range 


e J j 


1 (mK 2 ) 


a (/iK 2 ) 


x (/iK 2 ) 


350-550 


470 


2250 


92 


-345 


550-650 


608 


1982 


92 


-303 


650-730 


691 


1879 


89 


-267 


730-790 


763 


2180 


111 


-239 


790-850 


823 


2391 


115 


-255 


850-910 


884 


1824 


90 


-164 


910-970 


944 


1427 


70 


-104 


970-1030 


1003 


1111 


57 


-18 


1030-1090 


1062 


1043 


54 


17 


1090-1150 


1122 


1143 


57 


34 


1150-1210 


1182 


1067 


51 


75 


1210-1270 


1242 


808 


46 


119 


1270-1330 


1301 


693 


43 


151 


1330-1390 


1361 


778 


47 


193 


1390-1450 


1421 


746 


46 


218 


1450-1510 


1481 


604 


41 


211 


1510-1570 


1541 


517 


41 


229 


1570-1650 


1618 


435 


34 


261 


1650-1750 


1713 


363 


30 


242 


1750-1850 


1814 


344 


32 


264 


1850-1950 


1898 


227 


33 


170 


1950-2100 


2020 


217 


31 


203 


2100-2300 


2194 


162 


31 


244 


2300-2500 


2391 


159 


43 


357 


2500-3000 


2646 


105 


45 


560 



Note. — Band multipole range and weighted value £ e ff, decor- 
related band-powers qs, uncertainty as , and log- normal offset xb 
from the joint likelihood analysis of the 10 ACBAR fields. The pos- 
itive trend in the log-normal offsets with increasing £ is due to the 
increasing contribution of instrumental noise to the error budget; 
the log- normal offset approaches infinity in the Gaussian limit. The 
PMN radio point source and IRAS dust foreground templates have 
been projected out in this analysis. 
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Fig. 6.— Th e ACBAR band- powers plotted with those from WMAP5 (|Hinshaw et all |2008) and the 2003 flight of 
BOOMERANG Q ones et afl [2006). The three experiments show excellent agreement in the region where they overlap. 



sic six parameters. The results are consistent with the fiat case, but with the standard geometrical degeneracy 
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Fig. 7. — ACBAR results on the high-£ anisotropics. Left: The ACBAR band-powers above I = 1000 plotted against the best-fit 
ACBAR+WMAP5 model spectrum. The latest CBI results at 30 GHz and the previous ACBAR results are also shown. The ACBAR 
band-powers for I > 1950 are consistently below the reported high-£ CBI band-power. Right: The likelihood distribution for the ratio of the 
"excess" power, observed by CBI at 30 GHz and ACBAR at 150 GHz. The solid line assumes that dusty proto-galaxies do not contribute 
to the ACBA R b and-powers, while the dash-dot line includes the predicted contribution of approximately 22 /xK 2 at £ = 2600 from these 
sources (see i|6.2l l. With the inclusion of dusty galaxies, the likelihood peaks at a higher power ratio (~ 10). The excess power for each 
experiment is the difference between the measured and model band-powers for each experiment in a flat band with I > 1950. The likelihood 
for a given power ratio is found from Monte Carlo simulations of the band-powers and uncertainties. The vertical dashed line represents 
the expected ratio (4.3) for the excess being due to the SZE, while radio foregrounds would correspond to a ratio of > 10. If the excess 
power seen in CBI is caused by non-standard primordial processes, the ratio will be unity (blackbody), indicated by the dotted line. It 
is considerably more likely that the excess seen by CBI is caused by either the thermal SZE or radio foreground contamination than a 
primordial source. 

relating Qk and f^A expressed through 9 leading to a 
near-degenerate tail to Qk < 0. 

The ACBAR spectrum includes band-powers at I > 
2000 where the signal due to secondary CMB anisotropies 
associated with post-recombination nonlinear effects 
should become significant. In particular, the thermal 
Sunyaev-Zel'dovich effect and the contribution of unre- 
solved radio sources and dusty galaxies will ultimately 
dominate over the primary anisotropy damping tail; the 
only question is at what multipole crossover occurs. 
Without including such secondary effects, the parameters 
we derive from the primary anisotropy power spectrum 
could be biased. To account for this, we have added to 
the primary anisotropy power spectrum (1) the SZE tem- 
plate power spectrum Pf z used in K07 and ([Bond et al.l 
12003 iQoldstein et al.ll200l which was derived from cos- 
mological hydrodynamics simulations and (2) an unclus- 
tered point source template, as in eq.[3l Each template is 
scaled by an overall amplitude parameter, gsz and q src , 
which we assume have uniform prior measures with a 
range much larger than required by the ACBAR data. 
The white-noise form for 2?| rc given by eq. [3j is appropri- 
ate for the statistically-averaged power of a distribution 
of unclustered sources. The clustering of radio sources 
is not a large effect, but we do expect sub- mm sources 
associated with dusty galaxies at lower flux levels to be 
clustered. As mentioned in § 16.21 in spite of great strides 
in sub-mm observations in recent years, significant un- 
certainties remain in source fluxes and clustering at 150 
GHz. Theoretical models suggest both will be impor- 
tant for a complete treatment, but the approximation 
adopted here should be sufficient for the ACBAR data 
set. In the parameter tables below, we show results in- 
cluding these secondary templates. We find the basic 
parameter central values and uncertainties change little 
whether we marginalize over either of the two template 



amplitudes or set them both to zero. 

The parameter constraints are obtained using a Monte 
Carlo Markov Chain (MCMC) sampling of the multi- 
dimensional likelihood as a function of model param- 
eters. The pipeline is based on the public ly available 
CosmoMC 15 package ([Lewis fe Bridle! 12002]) . CMB an- 
gular power spectra and matter power spectra are com- 
puted using the CAMB code ([Lewis et al J 12000). As 
described in Section [71 we approximate the full non- 
Gaussian band-powe r likelihoods wit h an offset log- 
normal distribution ([Bond et alj l2000t ) . Our standard 
CosmoMC results include th e effects of weak gravita- 
tiona l lensing on the CMB (|Seliakl 119961 : iLewis et all 
120001) . Lensing effects in the temperature spectrum are 
expected to become significant at scales t > 1000, hence 
it is important to include this effect when interpreting 
the ACBAR results. The major effect of lensing is a 
scale-dependent smoothing of the angular power spec- 
trum which diminishes the peaks and valleys of the spec- 
trum. Inclusion of lensing in the model improves the fit 
to the data for all experiment combinations. 

The typical computation consists of eight separate 
chains, each having different initial, random parame- 
ter choices. The chains are run until the largest eigen- 
value of the Gclman-Rubin test is smaller than 0.01 af- 
ter accounting for burn-in. Uniform priors with very 
broad distributions are assumed for the basic parame- 
ters. The standard run also includes a weak prior on 
the Hubble constant (45 < H < 90 km s" 1 Mpc" 1 ) 
and on the age of the universe (> 10 Gyrs), but these 
have negligible effects. We also investigate the influ- 
ence of adding Large Scale Structure (LSS) data from 
the 2 degree Fiel d Galaxy Redshift Survey (2dFGRS) 
(ICole et al.l I2005D and the Sloan Digital Sky Survey 
CSDSS) ([Tegmark et alJl2006h . When including the LSS 

15 http:/ /cosmologist.info/cosmomc 
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data, we use only the band-powers for length scales larger 
than k ~ 0.1 h Mpc" 1 to avoid non- linear clustering and 
scale-dependent galaxy biasing effects. We marginalize 
over a parameter b 2 which describes the (linear) biasing 
of the galaxy-galaxy power spectrum for L+ galaxies rela- 
tive to the underlying mass density power spectrum. We 
adopt a Gaussian prior on b 2 centered around b g = 1.0 



with a very large width. We have also tried restricting 
the width to 8b, 
insensitive to this width 



g — 0.3, but the cosmic parameters are 



8.2. Base Parameter Results 

The results for the basic spatially flat tilted ACDM pa- 
rameters are presented in Table |4] The confidence lim- 
its are obtained by marginalizing the multi-dimensional 
likelihoods down to one dimension. The median value 
is obtained by finding the 50% integral of the resulting 
likelihood function while the lower and upper error lim- 
its are obtained by finding the 16% and 84% integrals, 
respectively. The CMBall data combination includes 
the ACBAR results presented here and other CMB data 
sets with published band-powers and windo w functions: 
the W MAP 5 year angular power spectra ;Nol ta et al.l 
(120081). and for co mpar ison the WMAP 3 year spec- 
tra lHinshaw et all (120061); the CBI extended mosaic re- 
sults dReadhead et al.l l2004aft and p olarization results 
(jReadhead et alJl2004bl : |Sievers et al. 12005ft combined in 
the manner descri bed inlSievers et alJ (12005ft : 16 the DASI 
two year results fHalvcrson et al. 2002ft: the DASI EE 
and TE band-powers (iLeitch et al.l FoOS); the VSA final 
results ([Dickinson et al.ll2004ft : the MAXIMA 1998 flight 
results (jHananv et al.ll2000ft : and the TT, TE, and EE 
resul t s from the BOOMERANG 2003 flight dJones et all 
120061: iPiacentini et al.ir2006t iMontrov et alJ [2006ft Only 
£ > 350 band-powers are included for BOOMERANG be- 
cause of overlap with WMAP (although inclusion of the 
lower £ results leaves the parameter results essentially un- 
changed). While ACBAR and BOOMERANG are both 
calibrated through WMAP, this is a small contribution 
to the total uncertainty in the ACBAR calibration and 
we treat the calibration uncertainties as independent in 
our parameter analysis. Although the DASI, CBI, and 
BOOMERANG 2003 EE and TE results for high-^ po- 
larization are included, they have little impact on the 
values of the parameters we obtain. 

The latest WMAP l ikelihood code found at 
http://lambda.gsfc.nasa.gov/ has been used in our 
analyses! When this ACBAR paper was submitted, 
all parameter analy ses were done using WMAP3 
(jHinshaw et all l2006ft : the new WMAP 5 results came 
out a few months later. The WMAP5 data allowed 
for an improved cross calibration of the ACBAR and 
WMAP band-powers and resolved the (< 1-er) parameter 
tensions that existed between the WMAP3 and ACBAR 
results. We note that the W MAP5 team's parameter 
analysis (jDunklev et al.ll2008l ) made use of the ACBAR 
band-powers with the 1.4% higher temperature calibra- 
tion from the WMAP3 and ACBAR cross-calibration. 
The marginalized one-dimensional likelihood distribu- 

16 We exclude the band-powers below £ = 600 from the CBI 
extended mosaic results to reduce the correlation with the TT 
band-powers of the CBI polarization dataset which influence the 
sample-dominated end of the spectrum. 



tions for the basic parameter set we obtain are shown 
in Figure [8] Note the contrast between the parameter 
determinations for WMAP 3 and WMAP 5. The primary 
improvement of WMAP5 over WMAP3 was a better 
understanding of the beam, which resulted in an im- 
proved measurement of the third acoustic peak. Other 
improvements in WMAP 5 were updated point-s ource 
correction, stimulated by (jHuffenberger et al]|2006ft . and 
foreground marginalization on large angular scales. 

The addition of LSS data has little impact on the mean 
values and errors in the cosmic parameters. The largest 
shift is < l-a in (l c h 2 , from 0.111±g;g^ to 0.107±^t- 
Our LSS results only include information on the shape 
of the density power spectrum, not its overall amplitude 
since we marginalize over the galaxy bias factor. If re- 
sults from weak lensing are included in the LSS data, 
then there is a slight increase in eg, but it depends some- 
what on which lensing results are included. The effect is 
to slightly increase the CMB+LSS result and improve 
the consistency with the CMBall results. The third peak 
is well determined by both ACBAR and WMAP5, and 
this defines the dark matter density; the inclusion of the 
LSS data does not significantly improve the constraints. 
We have also found that the results do not change signif- 
icantly from CMBall +LSS when SNla d ata are included 
(in this case from the lRiess et al.l (|2004ft gold set), so we 
have not included a separate column. Including SNla 
data would be crucial if we were attempting to constrain 
the equation of state of dark energy. 

All parameter results listed in Tables |4] and [5] include 
the effects of weak lensing of the CMB on the resulting 
power spectrum. For every case including CMB lensing, 
we have performed an identical calculation neglecting the 
effects of lensing. Including lensing improves the fit of the 
model to the observed band-powers compared for all data 
combinations. This can be quantified by the log-ratio 
of the lensed to no-lensed Bayesian evidence, Aln£ = 
In [_P( lens | data, theory) / ' P '(no — lens|data, theory)]. The 
evidence P (lens | data, theory) is an integral of the prod- 
uct of the a priori probability (the parameters' measure) 
and the likelihood of data given those parameters; it ap- 
pears in the denominator in the Bayesian chain to ensure 
the a posteriori probability has unity normalization. The 
resulting number is a conditional probability given the 
data and the assumptions about the parameters. The 
parameters and their measures are exactly the same, so 
the ratio is a robust indicator of preference. For WMAP5 
alone it is A In £ = 2.04; it increases to 2.89 with ACBAR 
included; and is 2.63 for CMBall. Naively relating this 
to a Gaussian translates to a significance of ~ 2.3cr for 
CMBall. From Fig. [5j it is clear that the difference in 
the power spectra of lensing and no-lensing is small; for 
this plot, the best-fit parameters from the lensed analysis 
were fixed and used to compute a lensed and unlensed 



spectrum. The difference, AC; 



lens 



c 



lens 



-lens 



IS 



explicitly shown in the inset of Figure We find only 
small shifts in the median value of the cosmic parameters 
when lensing is included; e.g., for the ACBAR+WMAP5 
data combination, we find er 8 = 0.79 



+0.03 
0.03 



0.80 



-0.03 
-0.03 



and fl c h 2 = 0.109i£:oo6 ~> O-H^aooe when g° in g from 
non-lensed to lensed models respectively. 

We now test whether the strength of the lensing mod- 
ification is consistent with expectations for lensing, by 
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Fig. 8. — Basic parameter marginalized 1-dimensional likelihood distributions for the following data combinations; WMAP3-only (green, 
dashed), WMAP5-only (black, solid), ACBAR + WMAP5 (red, long-dashed), CMBall (blue, dot-dashed). All runs include lensing. 



multiplying the lensing template AC] cns , which varies 
with cosmic parameters, by a variable strength q\ ens : 



^flens £j 



no— lens 



<7ic 



lens 



(5) 



The normalization is such that qi ens = 1 gives the nor- 
mal lensed CMB spectrum, while qi cns = gives the no- 
lensing case. An accurate determination of this subtle 
effect requires highly accurate window functions for the 
bands. We use a flat prior probability for gi en s, allowing 
it to vary from to 10. With WMAP5 alone, we obtain 



qic 

1.23 



1.34 



+0.21(+0.83). 
-0.23(-0.76)> 



+0.27(+1.51). 
-0.26(-0.85)' 



WMAP5+ACBAR gives 9lcns = 



CMBall gives q lcns = 1.21 



-0.24(- 
-0.24(- 



0.82) 
0.76) - 



We have also listed the 2-a errors (in brackets), which are 



far from twice the 1-sigma values, reflecting the highly 
non-Gaussian nature of the marginalized likelihoods ev- 
ident in the figure. Although we emphasize that qiens 
is not in any sense an independent parameter, it does 
illustrate that lensing of the expected strength is pre- 
ferred. 17 The significance of the detection is somewhat 

17 ICalabrcsc ct al. (2008) also undertook a lensing analysis of 
ACBAR temperature power spectrum, but, instead of qi cns as 
defined here, they used a multiplier Al of the lensing potential 
power spectrum, defined to be unity for normal lensing; they found 
A L = 3.0l° g for WMAP5+ACBAR. Repeating our analysis with 

this parameterization, we find lower values, Al = 1-60^q26(^o gg' ■ 
With the highest 3 ACBAR bins cut out, where secondary effects 
might have an impact on the determination of the lensing strength 
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Fig. 9. — 68% and 95% 2D marginalized contours in 6 and fi m for a number of data combinations. The results for 6 from WMAP are 
largely driven by the determination of the third peak. The value preferred by the WMAP5 data is slightly higher than that from WMAP3, 
and more consistent with that found when ACBAR is added to cither WMAP data set. 



less tha n the WMAP3/NV SS/S DSS cross-cor r elatio n re- 
sults of iSmith et all (|2007ft and iHirata et al.l (|2008l ). In 
Figure [TUl we show the marginalized distribution for the 
giens parameter using various combinations of CMB data. 

We have also run a limited set of non-flat model chains. 
The models in this case do not include the effect of weak 
lcnsing and we keep the same weak prior on Hq. When 
fife = is not enforced, the weak prior on Hq has a 
significant effect on the result as it restricts the extent of 
the geometrical degeneracy which is present in this case. 
For WMAP 5 only we obtain fl k = -0.018^;^ which 

becomes fi fe = — 0.013±g;g|| when ACBAR is included. 
8.3. Residual source marginalization 



As discussed in § 16.21 the ACBAR band-powers at 
£ > 1950 marginally exceed the predictions of the best-fit 
models for the primary CMB. We have repeated the ba- 
sic parameter runs including an unclustered source con- 
tribution described by equation [3] and marginalize over 
a wide uniform prior in g src from to 4600 fiK 2 , i.e. 
more than 100 times the power required to fit the high-i 
points. The results for runs including this marginaliza- 
tion are shown in Tables 2] and [5J The template ampli- 
tude <7 src is not well constrained by the fit, and the effect 
of the marginalization on the basic parameters is negli- 

parameter, the results are essentially unchanged. 



gible. The addition of a source template improves the fit 
to the high-^ points and therefore increases, albeit mod- 
estly, the best-fit likelihood: for the ACBAR+WMAP5 
combination the change in likelihood is AlnL = 0.24. 
The value and uncertainty for q src given in Table [4] spans 
the range of predictions for the contribution from dusty 
protogalaxies. 

Contributions from point sources or the SZE are nearly 
degenerate in the high-£ ACBAR band-powers. How- 
ever, the CMBall combination is potentially sensitive to 
an SZE contribution because of the SZE frequency de- 
pendence. The CBI data has a more significant excess at 
high-£, however, through substantial modeling and ve- 
toing, the CBI band-powers are expected to be free of 
radio source contamination. We assume that the resid- 
ual contribution of the unresolved source background to 
the CBI band-powers is sufficiently low that we do not 
require a source template for that data, only the SZE 
template. 

8.4. Sunyaev-ZeVdovich template extension 

The amplitude of the SZE signal depends strongly on 
the overall matter fluctuation amplitude, erg. We have 
modified our parameter fitting pipeline to allow for ex- 
tra frequency dependent contributions to the CMB power 
spectrum and have implemented it in a simple analysis 
using a fixed template 2?f z for the shape of the ther- 
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TABLE 4 
Basic 6 Parameter Constraints 
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Note. — Results for the basic parameter set. The runs all assume flat cosmologies, uniform and broad 
priors on each of the basic six parameters, and a weak prior on the Hubble constant (45 < Ho < 90 km s _1 
Mpc" 1 ) and the age of the universe (> 10 Gyr). Here A s = log[10 10 A s ]. All runs include the effect of weak 
gravitational lensing on the CMB. Column 5 presents the results when a SZ template characterized by the 
overall SZ-template-power gsz is included. The values of <rf z = q^Q (f2f,/i) 2//7 are higher than the erg derived 
from the primary anisotropics. Column 6 shows the results when the point source contribution scaled by 
q BTC is included in combination with the SZ-template with qsz set to al(flbh) 2 . The marginalization over 
either extra high-£ contribution does not significantly shift the results for the basic parameters. 



mal SZE power spectrum. The template was obtained 
from two large, hydrodynamical simulations of a scale- 
invariant (n s = 1) ACDM model wit h VL^h = 0.029 and 
with er 8 = 0.9 and a 8 = 1.0 . (See iBond et~aTl (pOOH ) 
for a detailed description of the simulations.) Recently 
the WMAP team have used a different SZE template 
based on analytic e stimations of the power spectrum 
(|Spergel et al.1 120061 ). It is characterized by a slower 
rise in t than the simulation-based template, which cut 
nearby clusters out of the power spectrum. There has 
been no fine-tuning of either spectra to agree with all 
of the X-ray and other cluster data. This may have an 
effect on shape, especially at high-£. 

The SZE contribution, X>f z = {q sz ) f v t>f z , added to 
the base six-parameter model spectrum has a frequency- 
dependent SZE pre- factor /„. Including this SZE tem- 
plate with all model parameters free to vary is comple- 
mentary to the analysis of § I7.2l which directly compared 
the residual CBI and ACBAR band-powers at i > 1950 
for the best fit WMAP5+ACBAR model power spec- 
trum. In that more restrictive analysis, the primary 
power spectrum is fixed and f u is allowed to vary as well 
as a broad-band excess power. We found the ratio of ex- 
cess power seen by CBI and ACBAR to be compatible 
with the ratio of the frequency prefactors /„ at 30 GHz 
and 150 GHz for the Sunyaev-Zel'dovich effect, although 
foreground contamination of both experiments could also 
contribute to the observed excess. In this section, we as- 
sume that all the excess power seen by the CBI experi- 
ment is due to the SZE. In the hydrodynamical simula- 
tions used to derive the SZE template, the amplitude was 
shown to scale as qsz = (a 8 /0.9) 7 (fl b h /0.029) 2 . We con- 
sider two cases: (1) the scaling parameter qsz is slaved to 
the results from the hydrodynamical simulations, which 



means that it is primarily determined by the primary 
CMB data; (2) qsz is allowed to float freely and an inde- 
pendent <7§ z is derived, to be compared with the as that 
is derived from the basic six parameters. We use a uni- 
form prior in qsz in this case with limits < qsz < 4.0. 

Regardless of the data combination, we find that in- 
cluding an SZE component in the model has little effect 
on the values of most basic cosmological parameters (see 
Table [4} , whether qsz is related to cosmic parameters 
through qsz = (cr 8 /0.9) 7 (Ob/i/0.029) 2 or is allowed to 
float freely. Note that the SZE results break the A s e~ 2T 
near-degeneracy (as does weak lensing, though not as 
strongly) . 

With the combination of the ACBAR and WMAP5 
data, for which there is only a weak indication of excess 
power, we find a freely-floating SZE amplitude results 
in qsz — 0.94lo;g3, with no effective lower bound. We 
can use the above relation of <7sz(cs, ^6^) to estimate a 
corresponding o'g SZ ' ) = 0.97^g "^g^Q g^ , where in brack- 
ets we have indicated the 2a errors. This result is higher 
than, but compatible within the uncertainties, to values 
obtained from the primary CMB fits alone: the ACBAR 
+ WMAP5 fits in Table SI give a s = 0.80^;^. The con- 

( sz) 

fidence limits of the derived <7g depend strongly on the 
choice of measure, which is here taken to be uniform in 
the amplitude qsz (This is evident in the translation of 
the relative flatness of the likelihood at low qsz to a "1- 



sigma detection" in a. 



(SZ), 



When the SZE contribution 



is slaved to the as and Vl^h values from the primary spec- 
trum, there is no effect on as- The high-^ excess power 
is not significant enough to change the well constrained 
value of us to the weakly preferred higher value. The 
effect of slaved and unslaved fits can be seen in Fig. [T3J 




Fig. 10. — The 1-dimensional marginalized likelihood distribution for the lensing amplitude qi cns for the combinations of CMB data 
shown. Although only qi ens = 1 is physically meaningful, these distributions indicate that the data prefer lensing of about the right 
magnitude. The ratio of the value at qi ens = 1 to that at qi cns = is the Bayesian evidence that lensing at the right amplitude is preferred 
over no-lensing. The inset figure shows AX>' cns = V\ cns - Dj, 10 ~ lcns , with both spectra computed using the best-fit parameters form the 
lensed analysis of the WMAP5+ACBAR data. Notice how small this correction is relative to D^ cns , shown as the dashed blue curve. 
The cosmological parameters used to construct AX>^ ens have been allowed to vary as the data requires in constructing the marginalized 
likelihoods for <?i cns . Each parameter variation leads to a slightly different difference template than that shown in the inset. 

When the high-£ band-powers of CBI and BIMA are 
included in the analysis, there is a significant detection 
of excess power. Both the CBI and BIMA band-powers 
are from 30 GHz interferometric observations and have 
higher /„ values than ACBAR. For the slaved case the 
as value and errors are unchanged. For the floating case, 

we find q sz = 0.691°;^ which maps to erf Z) = 0.93tg;^. 

These central values and uncertainties are computed 
by transforming integrals of the likelihood L(ggz, ^bh, ...) 



over the prior measures to confidence limits for (Tg 
Exactly what measure to place on q§i and therefore on 

a^^ is debatable. A measure uniform in a sz = q$z '■> 
as we used in K07 and (Goldste in et al.|[2003l ) translates 

— 1/2 

into a measure oc q sz dqsz which favors lower values 

of CTg than the measure uniform in qsz that we have 
adopted. 

As the cosmological parameters vary, the SZE tem- 
plate may depend on as and f^/i, and certainly depends 
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on the spectral index n s and astro-physical issues such 
as the history of energy injection into the cluster sys- 
tem. In a more complete treatment than that presented 
here, the shape should be modified along with the base 
cosmological parameters in the MCMC runs. 

We caution that the derived a g depends on the SZE 
template shape, its extension into the higher I regime 
probed by BIMA, and the prior measure placed upon 
gsz- 18 The modest decrease in log likelihood for the fit 
to the model when the SZE is taken into account, is 
Alni = 0.26 for the ACBAR+WMAP5 combination 
and increases to Ami = 1.91 for CMBall+BIMA. 

Regardless of the assumed prior, the interpretation of 
excess power as being due to the SZE results in a non- 
zero a g SZ%> for the CMBall+BIMA combination. Uncer- 
tainties for <7g are about a factor of two larger than 
for the as determined from the primary CMB data and 
there is a tension at about the 2-sigma level between the 
two median values. A visual summary of the results is 

shown in Fig. [T3l where we plot both as and a s against 
the spectral index for a number of data combinations. 
The addition of LSS data does not significantly change 
these results. When the ACBAR dusty point source or 
SZE template marginalization is included, as decreases 
slightly for the CMBall dataset. 

Recent weak lensing results are in basic agreement with 
the primary as values when J7 m = 0.26 ± 0.03 from Ta- 
ble H] is used. With 100 deg 2 of lensing data from the 
combination of the CFHT weak lensi ng legacy, RCS, 
Virmo s-Descart and GaBaDos surveys, iBeni amin et al] 
tf2007h get <7 8 (D m /0.26) a59 = °- 80 ± 0-05. From 
the C FHT weak lensing legacy survey alone, iFu et all 
120071) get <7 8 (f2 m /0.26) - 64 = 0.753 ± 0.043, and get 
o- 8 (fi m /0.26) - 53 = 0.82 ± 0.084 if only large-scale linear- 
regime results are used. These weak lensing numbers 
are lower than past published results because of im- 
proved treatments of the redshift distribution of the 
lensed sources. 

As shown in Fig. I12[ the marginal excess power in 
ACBAR is consistent with the combination of a point 
source contribution at the upper limit of the 150 GHz 
extrapolations and an SZE template at the level pre- 
dicted from the primary anisotropy as value. In this 
scenario, the cosmological parameters are virtually un- 
changed, but no explanation is provided for the CBI ex- 
cess. The SZE analyses with a free-floating amplitude 
have not included additional foreground sources for CBI, 
BIMA, or ACBAR. The effect of radio sources extrap- 
olated to 30 GHz was included in the original CBI and 
BIMA results, and is unlikely to be an important con- 
taminant for ACBAR. Dusty proto-galaxies should not 
effect CBI and BIMA, but will contribute to the high- 
£ ACBAR band-powers. Given the uncertainty in the 
source contamination for ACBAR, the weak detection 
of excess power does not significantly support the SZE 
interpretation of the CBI+BIMA excess. 

8.5. Running Spectral Index 

18 We also note that the non-Gaussian nature of the SZE signal 
was included in the BIMA results, but not in the CBI results. The 
non-Gaussian effect increases the sample varian ce and tends to 
open up the allowed range tow ards lower erg values (Goldstein et al. 
2003; Rcadhcad et al. 2004a). 



There has been interest in the running of the spectral 
index dn s /d\nk since the first release of WMAP data, 
which showed evidence for a significant negative running 
when comb ined with LSS and Lyman alpha forest ob- 
servations (jSpergel et al.l 120031 ). With the precision of 
the new ACBAR data, we might expect improved con- 
straints on running of the spectral index. To the basic 
six parameters in the minimal model, we add running of 
the spectral index dn s /d In k (fc*) around the pivot point 
fc* = 0.05 Mpc -1 . We adopt the conventionally-used 
uniform prior in dn s / d In fc(fc*), although in usual slow- 
roll-inflation models, the spectral index fluctuation 8n s cx 
ln(k/k i ,)dn s /dhik(k^) is typically restricted by |1 — n s \. 
Table [5] summarizes the parameter values when running 
is allowed and demonstrates that its inclusion has only 
a small effect on the other parameters. For WMAP5 
only, we find dn s /d\nk{k+) = -0.031 t° ° fs- The main 
tendency for the negative value comes from the low-£ 
WMAP5 data. When the ACBAR data is added the 
median value is similar, dn s /d\nk(k ir ) = — 0.037l Q22; 
with a reduced error. The results are nearly identi- 
cal for CMBall, but more compatible with no-running 
when LSS is added, — O.OlGj^ois- The precise measure- 
ment of the high-^ CMB power spectrum provided by 
ACBAR is a potentially powerful constraint on running, 
but may also include significant contributions from sec- 
ondary anisotropies and foregrounds. Therefore, in the 
last two columns we have included the effect of marginal- 
izing over the SZE or point source templates. In both 
cases, the mean value of the running becomes more neg- 
ative and more significant at the 2-sigma level. A vi- 
sual representation of the impact of adding the ACBAR 
data is given in Fig. [TT] which shows the correlation 
between n s and dn s /dlnk. The scalar spectral index, 
n s (K) = 0.918 ± 0.032 for the WMAP 5 + ACBAR 
combination, depends on the choice of pivot point k+; 
a smaller value would yield a higher result while a higher 
one would give an even lower result. 

8.6. Tensor modes 

We have run a limited number of cases including 
tensor modes, characterizing their strength relative to 
scalar perturbations by r and fixing the tensor tilt by 
n t ~ — r/8(l — r/16). The most stringent upper limit 
for r is given by the CMBall combination which yields 
r < 0.40 (95% confidence). This assumes a uniform prior 
measure for r, as is conventional in parameter estimation, 
although without much justification except that it is con- 
servative. Adding LSS tightens this lim it. Our result 
can b e compared with those obtained bv lDunklev et all 
(2008!), r < 0.43 (95% confidence) for WMAP5 alone. 

9. CONCLUSIONS 

We have used the complete ACBAR 150 GHz data 
set to measure the CMB temperature anisotropy an- 
gular power spectrum. Over three seasons of observa- 
tion, ACBAR dedicated 85K detector-hours to CMB 
observations at 150 GHz and covered 1.7% of the sky. 
The data are calibrated by comparing CMB tempera- 
ture maps for the largest ACBAR fields with those pro- 
duced by WMAP5. This calibration is found to be con- 
sistent with the previous planet-based and RCW38-based 
calibrations, but with the temperature uncertainty re- 
duced to 2.05%. In the original preprint of this paper, 
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TABLE 5 

Running Spectral Index Parameter Constraints 



WMAP5 



WMAP5+ACBAR 



CMBall 



CMBall+LSS CMBall+g sz CMBall+SZ+g sl 



n c h 2 



dn B /d ln(fc) 
As 
fi A 
Age 

Zre 

H 
qsz 

Qsrc 
°8 



o n2i q+ 0009 

U.UZlM_g OQOg 
n 117 +0.009 
u - 1± ' -0.009 
1 040+ - 003 

l.U^U_Q QQ3 

n 090 +a0()8 
u.uau_ 00g 

Q g2 o+0.043 

U.3ZO_ Q41 

_n n oi +0.029 
u.uoi_ 02g 
,+0.05 



3.09] 
0.7o] 
13.8 i 
0.30^ 



-0.05 
n +0.05 
-0.06 
= +0.2 
5 -0.2 
+0.06 
-0.05 

o.8il : ot 
n.7±i:J 

68.4^; 2 8 



+0.0007 
0.0007 
+0.008 
0.008 
+0.003 
0.003 
+0.009 
0.008 
+0.033 
0.031 
n o 7 +0.023 
U.VOI _Q.022 

3.10+ 005 



0.0219 
0.119 
1.041 
0.091 
0.918 



0.69 
13.8 
0.31 
0.83 
11.9 
67.9 



0.04 
+ 0.04 
-0.05 
+0.1 
0.1 
+0.05 
0.04 
+0.03 
0.03 
-1.7 
-1.6 
+3.5 
3.3 



022n+ 0006 
u.uzzu_ 0006 

120+ ' 008 

U.1ZU_Q QQ g 

i.u<iz_ 0Q2 
u.uyz_ 00g 
0.914lg;°i° 
-0.0391°;^ 

n ln +0 05 
°- lu -0.04 

0.691 
13.8j 
0.3i: 
0.83 H 



+0.0006 
-0.0006 
+0.005 
-0.005 
+0.002 
0.002 
+ 0.009 
-0.008 
Q48+0-026 

-o 016 +0019 

U.U1D_q 01g 

3.07+ 04 



0.0225 
0.109 
1.042 
0.091 



0.0217 



+0.0006 
0.0006 
+0.008 



n+0.04 
-0.05 
d+0.1 
-0.1 
,+0.05 
-0.04 
O+0.03 
-0.03 

i2.oti:S 

v+3.3 
-3.2 



67.7" 1 



0.75 
13.7 
0.25 
0.79 
11.5 
72.4 



0.04 
+0.02 
-0.03 
+0.1 
0.1 
+0.03 
0.02 
+0.03 
-0.03 
+ 1.6 
1.6 
+2.4 
2.4 



0.122_ 00g 

1 041+0-003 

± - u ^ -0.002 

092+ - 009 
u.uy^_ 009 

894+ - 031 

-0 052+ 022 
u.uoz_ 021 

-,+0.05 
-0.04 
j+0.04 

-0.05 

D + 0.1 

5 -0.1 
7 +0.05 
-0.04 
,+0.03 
-0.03 
,+1.8 
-1.7 
7+3.3 
-3.1 
7+0.12 
-0.11 



3.10] 
O^ 
13.! 
0.32]| 
0.83^ 
12.2^ 
66.7j 
0.77j 



0.93] 



021 7 +0 - 0007 
u - uzl ' -0.0006 

, 24+0.008 
u - iz ^-0.008 

1 041+0-002 

Lml -0 .002 

094+ 009 

n SQfi+0-030 
0-89b_ 030 

-0 052+ ' 022 
u.uoz_ 022 

, +0.05 

-0.04 

7+0.04 

-0.05 

D + 0.1 

-0.1 
n oo+0.05 
U ■ -3-3 _ o.04 
1+0.03 
-0.03 

12.41J-S 
-,+3:2 
-3.2 



3.11] 
0.671 
13.8]| 
0.331 
0.841 
12.4] 
66. 2 4 



30] 



Note. — Marginalized parameter constraints for models with a running of the spectral index. The tendency 
for the low-^ data to prefer negative values continues with the higher £ data, but at less than 2-sigma for CMB- 
only. The last two columns include marginalization over extra high-i? contributions from SZ and point sources. 
The basic parameters are stable with respect to this marginalization while the median value for dn s /dlnfe(fc*) 
is increased slightly when the SZ or point source contribution is included. 
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n. 



Fig. 11. — 68% and 95% 2D marginalized contours in n s and n rU n = dn s /d In k. The left panel shows the results for the WMAP5 and 
ACBAR+WMAP5 combinations. The right panel shows the effect of adding the rest of the CMB data and marginalizing over the SZE 
template. The basic parameters are essentially unchanged if we marginalize over one or the other of the two template amplitudes. 



the ACBAR band-powers were calibrated through com- 
parison with WMAP3. The new WMAP5 results be- 



came available during the review process and we have 
updated the ACBAR band-powers and cosmological pa- 
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Fig. 12. — Best-fit models for the ACBAR + WMAP5 combination. Only the I > 1000 ACBAR band-powers are shown. The arrows 
are indicative of the possible (coherent) shift in the 1-<t confidence limits due to beam error. The top panel shows fits with just primary 
CMB (black, solid) and with an extra SZE contribution. The case where the SZE template is scaled independently as gsz gives the best 
fit to the high-f excess (red, short-dashed) with ggz = 0.82. The case where the SZE amplitude is slaved to the cosmological parameters 
<7g(r2(,h) 2 (blue, long-dashed) does not yield enough power to fit the excess given the best-fit value of ag = 0.81. The bottom panel shows 
the case when a point source contribution scaled by q B ic is included for ACBAR together with the slaved SZE contribution (blue, dotted). 
In this case the best-fit model has ag = 0.81 determining the sub-dominant SZE contribution and q src = 31[iK 2 determining the point 
source contribution. 



rameters to take advantage of this new information. The 
ACBAR band-powers are otherwise unchanged. The 
new WMAP5 parameters and calibration subtly changed 
the best fit WMAP5 and ACBAR model resulting in a 
decrease of the significance of the high-£ excess in the 
ACBAR data. 

The ACBAR band-powers reported in Table [3] are 
the most sensitive measurements to date of CMB tem- 
perature anisotropies for multipoles between £ ~ 900 
and 3000. In this data, the fourth and fifth acous- 
tic peaks are significantly detected for the first time. 
These precise measurements of the CMB temperature 
anisotropies at high-€ are consistent with a spatially 
flat, dark energy-dominated ACDM cosmology. Includ- 
ing the effects of CMB weak lensing in the computation 
of model power spectra improves the fits to the combined 
ACBAR+WMAP5 data. The excellent fit of the ACDM 
cosmological model to the combined ACBAR+WMAP5 
data at I < 2000 is a strong confirmation of the stan- 
dard cosmological paradigm and gives us confidence in 



the resulting parameter values. The ACBAR data favor 
higher median values of as, f2 TO > and than those pre- 
ferred by WMAP3; however, with the improved WMAP5 
results for I > 650, this tension has been resolved. The 
parameter values remain stable with the inclusion of ad- 
ditional CMB and LSS data sets, with or without the 
marginalization over an SZE or point source template 
for ACBAR. For example, as ~ 0.80 holds for all param- 
eter variations in Table [U and even in Table with the 
inclusion of a running spectral index. 

We have performed strict jackknife tests with the data, 
and find that the results are free of significant systematic 
errors. We have projected out templates derived from 
the FDS99 dust model and the PMN radio source cata- 
log from the ACBAR maps before estimating the band- 
powers and find the residual contributions from these 
foregrounds to be negligible at the current sensitivity. 
The contribution of dusty proto-galaxies is expected to 
be insignificant for all but the few highest-^ band-powers, 
but remains poorly constrained due to our incomplete 
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Fig. 13. — The figure contrasts the one and two sigma contour intervals for erg determined from the primary anisotropy component of 
the CMB (left) with the value inferred from the SZE template transformation of <jgz into <Xg (right), assuming a uniform prior measure 
in qsz- Allowing for a point source contribution would decrease the tension between erg and a^ z for the ACBAR+WMAP5 case. These 
panels also demonstrate the strength of the deviation of n a from unity for the fiat ACDM model. 



knowledge of the spectral dependence of these sources. 

Secondary anisotropies are expected to become impor- 
tant at small angular scales. The ACBAR band-powers 
are slightly higher (1.1a) than expected for the primary 
CMB anisotropy at multipoles above t ~ 2000. We ex- 
pect that some of this power results from contamina- 
tion by dusty proto-galaxies. However, the combined 
signal is considerably smaller than the significant detec- 
tion of excess power reported by the CBI experiment at 
30 GHz. A joint analysis of the CBI and ACBAR band- 
powers in the multipole range of 2000 < I < 3000 argues 
strongly against the CBI excess having the spectrum of 
primary CMB anisotropy. These results are consistent 
with the excess power seen by CBI being due to either 
the Sunyaev-Zel'dovich effect or radio source contamina- 
tion. Higher sensitivity observations over a broad range 
of frequencies are necessary in order to fully character- 
ize CMB secondary anisotropies and eliminate potential 
foreground contamination. 
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TABLE 6 

Error Budget for the a£ m -BASED ACBAR Calibration 



Source 



Uncertainty (%) 



Statistical Error in the Calibration ratio 



1.10 
1.1 
0.35 
0.5 
1.0 
0.3 
0.5 
0.2 



t dependence of the Calibration ratio 
Statistical Error in the Transfer Function of the Calibration ratio 
Uncertainty in the WMAP B e 
Relative pointing uncertainty 

Uncertainty in the Year-to-year ACBAR calibration 
Uncertainty in the Transfer Function for the Power Spectrum 
Contamination from foregrounds 



WMAP5's Absolute Calibration 
Overall 



0.2 
2.05% 



Note. — The calibration of ACBAR using the WMAP5 temperature maps 
has multiple potential sources of error, tabulated here for reference. The domi- 
nant calibration uncertainties are due to noise in the WMAP maps at the angu- 
lar scales used for calibration. The uncertainty in the ACBAR beam function 
is comparable to the calibration uncertainty. 



APPENDIX 



CALIBRATION 



The calibration used in K07 was linked to the Boomerang03 (B03) calibration with observations of RCW38. In 
this section, we describe a new calibration using an a/ m -based comparison of CMB structure observed by WMAP5 
and ACBAR in 2005. This method is inspired by the calibration scheme used to calibrate B03 to WMAP. The 2005 
calibration is carried to other years by an ACBAR-ACBAR power spectrum comparison on fields observed in both 
years. The WMAP- ACBAR cross-calibration method is described below, with a detailed accounting of uncertainty in 
Table [6] 

WMAP-ACBAR Calibration — Calibrating with the CMB temperature anisotropies has two main advantages. The 
first is that the calibration of the WMAP temperature maps (at 0.2% in temperature) is an order of magnitude more 
precise than the flux calibration of the calibration sources ACBAR used in previous releases. The second advantage is 
that the anisotropies have the same spectrum as what is being calibrated, rendering the large frequency gap between 
WMAP and ACBAR irrelevant. 

The two experiments have different scan patterns, noise, beam widths, and spatial filters that will effect the measured 
flux. In this analysis, we assume that the WMAP5 maps are effectively unfiltered except for the instrumental beam 
function. The two maps can then be represented as: 



where T is the underlying CMB signal, N is the instrumental noise, B is the beam function, and Fij is the ACBAR filter 
matrix as defined in Section [31 We reduce the filtering differences by resampling the WMAP map using the ACBAR 
pointing information and applying the ACBAR spatial filtering to generate an 'ACBAR-filtered' WMAP map. 



The results of applying this algorithm to the B03 map is shown in Figure [TU We choose to do the absolute calibration 
via cross-power spectra rather than a direct pixel-to-pixel comparison of the maps. Using cross-spectra significantly 
reduces the impact of the noise model on the result. The significant beam differences between the experiments are 
more naturally dealt with in multipole space than in pixel space. We construct the ratio from the filtered maps: 



where X can denote either the V- or W-band map for WMAP and Y/Z marks either of two noise-independent ACBAR 
combinations. There is a narrow Grange from 256-512 useful for calibration. The range is limited at high-£ by the 
rapidly falling WMAP beam function and at low-£ by the ACBAR polynomial filtering which acts as a high-pass filter. 
We choose to use the WMAP V & W bands to take advantage of their smaller beam size. 

Monte Carlo simulations are used to determine the transfer function of this estimator. We generate CMB sky 
simulations convolved with the respective instrumental beam functions using the Healpix 19 library. We resample 

19 http://healpix.jpl.nasa.gov 
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Fig. 14. — A comparison of observations of the CMB8 field made with B03 and ACBAR. This field lies in the deep region of the B03 
map. The top two maps are from B03. The bottom two maps are from ACBAR. In the top left panel, the B03 map of the CMB8 field. The 
dynamic range of this map is greater than that of the other three figures. The increased noise at one edge marks the edge of the B03 deep 
coverage. The ACBAR filtering is applied to the B03 map to create the map in the top right panel. Directly below it in the bottom right 
panel is the ACBAR map of same region. Note the clear correspondence between the CMB anisotropics observed by B03 and ACBAR. 
Three bright point sources have been masked. An ACBAR left-right sweep difference map is shown in the bottom left panel. The power 
spectrum of this map (and the other 9 fields) is plotted in Fig. \3\ 

each realization and apply the ACBAR filtering matrix described above to generate equivalent maps for each field. 
We expected and found a small intrinsic bias as the beam convolution and filtering operations do not commute: 
B AC bar t f .. b wmap ^ B WMAP * F .. B ACBAR We correct the rea i da ta by the ^-dependent transfer function 

measured in these simulations. The technique is easily adapted to estimate the error caused by pointing uncertainties 
and to confirm that the estimator is unbiased with the inclusion of noise. The derived error in the transfer function is 
listed in Table H 

Foreground sources have the potential to systematically bias a calibration bridging 60 to 150 GHz. Radio sources, 
synchrotron emission, dust, and free- free emission all have a distinctly different spectral dependence than the CMB, 
which could lead to a calibration bias. This risk is ameliorated by the positioning of the ACBAR fields in regions of 
exceptionally low foregrounds. Bright radio sources detected in either experiment are masked and excluded from the 
calibratio n. The calibration pro ved insensitive to the exact threshold for source masking. We use the MEM foreground 
models in Hin shaw et al.l (|2006f ) to estimate the RMS fluctuations of each foreground relative to the CMB fluctuations 
and find that the free-free and synchrotron fluctuations are less than 0.1% of the CMB fluctuations in all frequency 
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bands while dust emission can reach 1.5% of CMB fluctuations in the 150 GHz maps. We test the effects of the most 
significant foreground, dust, by adding the FDS99 dust model (jFinkbeiner et al1l!999T l to a set of CMB realizations. 
The resultant maps arc passed through a simulated pipeline as outlined in the previous paragraph. We find that the 
addition of dust does not introduce a detectable bias with an uncertainty of 0.2%. 

We perform a weighted average of the measured calibration ratio across all ^-bin, field and band combinations after 
correcting for the estimated signal-only transfer functions. We estimate the calibration error to be 1.97% for the 2005 
data. Table [S] tabulates the contributing factors and error budget. We then propagate this a^ m -based calibration to 
the CMB observations done in 2001 and 2002. 

ACBAR 2001-2002 and 2002-2005 Cross Calibrations— We propagate the 2005 calibration into 2001 and 2002 by com- 
paring the 2001 observations of the CMB2 field to the overlapping 2002 CMB4 field, and the 2002 observations of the 
CMB5 field to the 2005 observations of the CMB5 field. A power spectrum is calculated for each overlapping region 
and the ratio of the band-powers is used to derive a cross calibration. The procedures used are outlined in more detail 
in K07. We use the same relative calibration for 2001 as K07: T2001/T2002 = 1.238 ± 0.067. We find cross-calibration 
factor for 2002 to be T2005/T2002 = 1.035 ± 0.025. We apply these corrections to the data and determine the over- 
all calibration uncertainty to be 2.05% (in temperature units) based primarily on the uncertainties associated with 
WMAP/ACBAR-2005 cross calibration. 
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